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Introduction and Summary of Paper. 

Problems which deal with the stabihty of bodies in equilibrium under stress are so 
distinct from the ordinary applications of the theory of elasticity that it is legitimate 
to regard them as forming a special branch of the subject. In every other case we 
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are concerned with the integration of certain differential equations, fundamentally 
the same for all problems, and the satisfaction of certain boundary conditions ; and 
by a theorem due to Kirchhoff"^ we are entitled to assume that any solution which we 
may discover is unique. In these problems we are confronted with the possibility of 
two or more configurations of equilibrium, and we have to determine the conditions 
which must be satisfied in order that the equilibrium of any given configuration may 
be stable. 

The development of both branches has proceeded upon similar lines. That is to 
say, the earliest discussions were concerned with the solution of isolated examples 
rather than with the formulation of general ideas. In the case of elastic stability, a 
comprehensive theory was not propounded until the problem of the straight strut had 
been investigated by Eulek/j" that of the circular ring under radial pressure by 
M. Levy J and G. IT. Halphbn,§ and A. G. Green"HILL had discussed the stability of 
a straight rod in equilibrium under its own weight, || under twisting couples, and when 
rotating. II 

In a paper which has become the foundation of the theory in its existing form,'^''^ 
G. H. Bryan has brought these isolated problems for the first time within the range 
of a single generalization. Examining the conditions under which Kirchhoff's 
theorem of determinacy may fail, he was led to the conclusion that instability is only 
possible in the case of such bodies as thin rods, plates, or shells, and in these only 
when types of distortion can occur which do not involve extension of the central line 
or middle surface, so that it is legitimate to discuss any problem in elastic stability 
by methods which have been devised for the approximate treatment of such 
bodies. He showed, moreover, that the stability of the equilibrium of any given 
configuration depends upon the condition that the potential energy shall be a 
minimum, in that configuration. 

A closer examination of Bryan's theory suggests that some of the conclusions 
which have been drawn from it are scarcely warranted. The contention that no 
closed shell can fail by instability, because any distortion would involve extension of 
the middle surface, will be discussed later. ft For our present purpose it is sufficient 
to remark that the whole theory is based upon the assumption that the strains 
occurring previously to collapse must be kept to the extremely narrow limits within 
which, in the case of ordinary materials, Hooke's Law is satisfied. This assumption, 
of cou^rse, expresses a restriction necessarily imposed upon the range of practical 

'^ A. E. H. Love, 'Mathematical Theory of Elasticity' (second edition), § 118. 
t 'Hist. Acad. Berlin,' XIIL (1757), p. 252. 

I 'Liouville's Journal,' X. (1884), p. 5. 

§ 'Comptes Rendus,' XCYIIL (1884), p. 422. 

II 'Proc. Camb. Phil. Soc.,' IV. (1881), p. 65. 
^ 'Proc. Inst. Mech. Eng.,' 1883, p. 182. 

•^^ 'Proc. Camb. Phil. Soc.,' VI. (1888), p. 199. 
tt Cf, pp. 222, 236, 
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problems which can be treated by the ordinary theory of elasticity ; but it is not 
legitimate to conclude that instability is only possible, even if its conditions were only 
calculable, in the case of materials which obey Hooke's Law, and there is no warrant 
for the employment of '^ crushing formulae" in the design of short struts and thick 
boiler flues/^' 

A more serious weakness in the existing theory of elastic stability, when regarded 
from the mathematical standpoint, is the fact that the methods which it employs are 
admittedly only approximate. The higher the elastic limitf of the material under 
consideration, the less adequate are these methods to deal with the whole range of 
problems which should come within its scope. In fact, we are faced with the 
anomaly that, while in its ordinary applications the theory of elasticity is not 
concerned with the conception of an elastic limit, in questions of stability the 
existence of finite limits is an essential condition for the adequacy of its results. In 
an ideal material, possessing perfect elasticity combined with unlimited strength, 
types of instability could occur with which existing methods would be quite 
insufficient to deal. 

The theory of elastic stability is thus in much the same position as that of the 
ordinary theory of elasticity before the discovery of the general equations, and one 
aim of the present paper is to remedy its defects by the investigation of general 
equations, which may be termed '' Equations of Neutral Equilibrium," and which 
express the condition that a given configuration may be one of limiting equilibrium. 
These equations are universally applicable only to ideal material of indefinite strength, 
and the possibility of elastic break-down must receive separate investigation ; but 
they are also applicable, even with materials of finite strength, to any problem which 
comes within the restrictions imposed by Bryan's discussion, and therefore enable us 
to test the accuracy of his treatment of problems, such as that of the boiler flue, for 
which the ordinary Theory of Thin Shells has been thought insufficiently rigorous.J 

In every problem of this paper it is found that the Theory of Thin Shells gives a 
solution which is correct as a first approximation, and the practical advantages of the 
new method of investigation are, therefore, not immediately apparent. But it must 
be remembered that the approximate theory of thin plates and shells has not as yet 
been rigorously established, and that much work has recently been undertaken with 
the object of testing it by comparison with accurate solutions of isolated problems.§ 
Now in finding conditions for the neutrality of the equilibrium of any given 
configuration we are at the same time obtaining the solution of a statical problem ; 
for a configuration of slight distortion from the equilibrium position will also be one 

^ W. C. Unwin, 'Elements of Machine Design' (1909), Part L, p. Ill ; S. E. Slocum, "The Collapse 
of Tubes under External Pressure," 'Engineering,' January 8, 1909. 

t By " elastic limit " is intended, here and throughout this paper, the limit of linear elasticity. 

J Cf. pp. 210, 224. 

§ Love, op, cit., Introduction, p. 29, and Chapter XXII. 
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of eqiiilibrmm. Hence every solution which we can obtain will add to the number 
of these ^' test cases," which has not hitherto included solutions for any but plane 
plates. 

A far more important advantage of the new method , from the practical point of 
view, is the accuracy with which it follows the actual "stress history" in a body 
which fails by instability under a gradually increasing stress. In cases where 
instability precedes elastic break-down this difference of method is not important ; 
but for the discussion of instability in overstrained material, where the stress-strain 
relations are intimately dependent upon the previous stress history, its introduction 
is absolutely necessary. 

The extension of Euler's theory to struts of practical dimensions and materials, 
which forms the conclusion of this paper, suggests a large and new field for 
investigation. The number of similar cases which can be treated, in the existing 
state of our knowledge of plastic strain, is very small, and indications are given below 
of the questions which still require an answer ; there is reason to believe that the 
requisite experimental research would not present insuperable difficulties, and that 
we may hope in the future to obtain an adequate theory of experimental results 
which are at present very little understood. 

Equatio:n's of Netjtkal Equilibrium in Rectangular Co-ordinates. 

Method of Derivation. 

The question of stability arises in regard to any system in which there is a 
possibility of slight displacement from the configuration of equilibrium. This possi- 
bility may be afforded either by a more or less limited degree of mechanical freedom— 
in which case the problem is one of statical stability, and practically unaffected by 




Fig. 1. 

the tendency, which any actual body displays, to distort under the infiiuence of 
applied forces ; or it may be due, more or less entirely, to this tendency. In the 
latter case the problem is one of elastic stability, and must be treated by distinct 
methods. There is, however, no essential difference between the two types of 
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instability, and a general discussion of the elastic type may be very conveniently 
illustrated by reference to a mechanical example. 

In this connection we may consider the system illustrated by fig. 1, in which a 
uniform heavy sphere rests in equilibrium within a hemispherical bowl, under the 
action of its own weight and of the pressure exerted by a pointed plunger, which is 
free to move in a vertical line through the centre of the bowl. This system has been 
chosen for the illustration which it affords of collapse under a definite '^ critical 
loading." In this it bears an unusual resemblance to examples of elastic instability — 
the stability of most mechanical systems being dependent solely upon the relative 
dimensions of their members. In the absence of friction, we find that the equilibrium 
will become unstable as the load on the plunger is increased through a critical value 
given by -yrr 

W is the weight of the sphere, 
r is the radius of the sphere, 
R is the radius of the bowl. 

The above solution rests upon the assumption that the sphere, bowl and plunger 
are absolutely smooth and rigid, and the possibility of slight displacement is afforded 
by the freedom of the sphere to take up any position of contact with the bowl. To 
discuss the equilibrium of the sphere in the position illustrated we must consider the 
forces which act upon it in a position of slight displacement. These include two 
systems, one tending to restore the initial conditions, the other tending to increase 
the distortion, and stability depends upon the relative magnitude of the two effects. 
We may investigate the problem by three methods, fundamentally equivalent, which 
are described below : — 

(1) The Energy Method. — We may derive expressions for the potential energy of 

the system in a position of slight displacement from the equilibrium position. 
The condition of stability requires that the expression for the potential energy 
shall have a minimum value in the equilibrium, position. 

(2) The Method of Vibrations. — We assume that the slight displacement has been 

effected by any cause, and investigate the types of vibration possible to the 
system when this cause is removed. The condition of stability requires that 
all such types shall have real periods. 

(3) The Statical Method. — We confine our attention to the special case in Avhich 

the stability of the equilibrium position is neutral. In this case there must 
exist some type of displacement for which the collapsing and restoring 
effects, discussed above, are exactly balanced, so that it may be maintained 
by the original system of applied forces. We have, therefore, to find 
conditions for the equilibrium of a configuration of small displacement, under 
the given system of applied forces. 
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Any of these methods is vaHd for the investigation of elastic stability, and all have 
in fact been employed, the displacement considered being that of the central-Hne or 
middle-surface of the rod or shell, and the resultant actions over cross-sections being 
derived in terms of this displacement, by the approximate theory first suggested by 
KniOHHOFF. The third method is generally found to be preferable, and is the basis 
of the investigation to be described below, but the actual procedure will be found to 
possess one or two novel features. 

In 'the first place, an endeavour will be made to dispense with the assumption that 
elastic break-down occurs at very small values of the strains ; instead, we shall deal 
with an ideal material possessing perfect elasticity combined with unlimited strength. 
Such a material could not faily unless hy instability^ and our problems will no longer 
be confined to thin rods, plates, or shells. It follows that we can only obtain sufficient 
accuracy in our conditions for neutral stability by deriving them with reference to a 
volume-element of the material. 

Further, since instability will in some cases not occur until the strains in the 
material have reached finite values, we shall have to introduce an unusual precision 
into our ideas of stress and strain. The discussion of finite strain is merely a problem 
in kinematics, and has been worked out with some completeness"^ ; but the corre- 
sponding stress-strain relations in our ideal material are necessarily less certain, since 
they must be based upon experiments in which only small strains are permissible. 

For example, if we assume that Hooke's Law is satisfied at all stresses, we must 
decide whether our definition of stress is to be 



Lt. 

or 

Lit, 



Total action over an element of surface 
Original area of that surface 

Total action over the surface 



„Area of that surface after distortion 

For the ordinary purposes of elastic theory the two definitions may be regarded as 
equivalent, and the distinction is too fine to be settled experimentally. In the 
absence of any generally-accepted molecular theory which might indicate the correct 
result, it seems legitimate to make the simplest possible assumptions which do not 
involve self-contradictions, and which yield the usual results when the strains are 
very small. 

It may be shown t that in a distortion of any magnitude three orthogonal linear 
elements issue from any point after distortion, which were also orthogonal in the 
unstrained configuration, and that these linear elements undergo stationary {maximum 
or minimum or minimax) extension. Hence an elementary parallelepiped constructed 
at the point, with sides parallel to these linear elements, undergoes no change of 
angle in the distortion. It is clear that only normal stresses will act upon its faces 

"^ For a discussion of the theory, with references, see Love, op. cU., Appendix to Chapter I. 
t Love, op. ciL, §§ 26, 27. 



MR. R. V. SOUTHWELL ON THE GENERAL THEORY OF ELASTIC STABILITY. 193 

after distortion, and that if these stresses be expressed in terms of the extensions of 
the sides we have complete relations between stress and strain. 

We shall therefore assume that these principal stresses and principal strains, 
whatever their magnitude, are connected by the ordinary equations of Hooke's Law ; 
that is to say, if the extensions in the principal directions are Ci, 62, e^, and the 
corresponding stresses are Rj, R^* R;]> then 

1 



e, = 



m 



» * • • » CX/v^. > 



1 

where E is Young's Modulus, and — is Poisson's ratio for the material under 

m 

consideration. 

These relations may be written in the form 

^1 = 7 — -TYT ^[(w.-l)ei + e2 + ^J 

(m+l)(m — 2) 

[(m— 1) 61 + ^2 + ej, ..., &c., (2) 



m — 2 

where G is the Modulus of Rigidity. 

In these relations the measure of extension is assumed to be 

Increase in length of linear e lement 

Length of the element before strain ' 
and of stress^ 

Total action over an element o f surface 
Area of the element before strain 

We have then the usual expression*!* for the energy of strain, per unit volume of the 
unstrained material, in terms of the principal extensions, viz. : — 

W ==^^C! (61 + 62 + 63)^-^20 {6162 + 62^3 + 63^1). ..... (3) 

III "~" Zj 

The above assumptions yield sufficient data for the calculation of the stress system 
in any configuration of equilibrium, even when the strains are not small. Assuming 
that the calculation has been effected, we have to show how conditions for the 
stability of the system may be obtained. 

We must distinguish three configurations : the unstrained configuration, in which 
the co-ordinates of any point are given by x, y, z; the configuration of equilibrium 
under the stress- system, the stability of which we are investigating ; and a configu- 

■^ This assumption is open to the objection that it would render possible the compression of a material 
to zero volume by means of a finite stress. It will nob, however, introduce any serious error, and has the 
advantage, which more probable assumptions do not possess, of leading to a definite energy-function. 
The definitions of stress and strain given above are generally employed in the construction of " stress- 
strain diagrams " from a tension test, the extensions of the specimen being taken as abscissae, and the 
total loads as ordinates of the plotted curve. 

t Love, ojx ciL, § 68. 

VOL. CCXIII. — A. 2 C 
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ration of slight distortion from the equiKbrium position, which can be maintained 
vyithout the introduction of additional stress at the boundaries, if the equilibrium of 
the second configuration is neutral. We shall consider first a stress-system which is 
such that the principal stresses in the second configuration have the same magnitudes 
and directions throughout the body ;^ and we shall take these directions as axes of 
X, y and z. We may then define the second and third configurations by saying that 
in them the co-ordinates of the point (^r, y^ z) become 



an 



d 



Cc(l+(:^l), 



^/(l+eg), 



z 



(1+63), 



^ ( 1 + ^1) + ^^ 1/(1 + 62) + v\ ^ { 1 4- eg) + 1//, 



respectively. We shall not limit the values of ej, e^, e^, although in practical cases 
they must be small : u' ^ v\ w' are infinitesimal. In the second configuration the axes 
Ox^ Oy^ Oz are directions of principal stress, and the stresses are 



jL 



X 



20 

m — 2 



[(m — l)(3i + 62 + ^3], ..., &c. . ..... {2) his 



In the third configuration we shall find that lines which in the first configuration 
were slightly inclined to Ox, Oy, Oz become directions of principal stress and strain. 
The final extension of a line which originally had direction-cosines l, m, n is 



e 



r "T* 



V 



lll-^e,+ ^^'^-^ 



did , 3t/^ 
+ in -r h n 



2 



ox) ' "^ oy ' " dz 



+ U- — hm 1 +60+ ;-- +r^ -;-" r + 1 ^ "^ r-m -7^ — hn i+e.^+ -tt- 

t dx \ oy I oz] L dx oy \ oz 



. . . (4) 



It may be shown that e! has a stationary value when 



A 



m ~ nil = 



an 



d 



n = Ui 



(l+eO'-(l+e,)^ ' 

dz dx_ 



» . . \ / 



J 



to terms of the first order in %i\ v\ %v\ 

^ In some cases, sach as Greenhill's problem of tlie stability of a heavy vertical rod (p. 188, footnote), 
it is necessary to allow for variation in one or more of the principal stresses ; the necessary alterations are 
easily made, and as they are nob required for the examples of this paper their consideration would involve 
unnecessary complexity. 

t Added May i.— The approximation of these expressions is insufficient if any two of the principal 
strains (ei, 62, e^) in the second configuration are equal ; in this case additional terms must be retained in 
the denominators. The equilibrium under hydrostatic stress (^<i = ^2 = e-s) is necessarily and obviously 
stable.] 
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* 

Thus the line initially given by the direction-cosines 

becomes a direction of principal stress in the final configuration. Its direction-cosines 
(referred to Ox, Oy, and Oz) are then 



or 



1 



1, 



at/ 

ax 



+m^ (1+^2) 



1 -f- e. 



c)x 



{l+e,y-{l+e,y 



CJiAy 



-{-fhi (14-63) 



1 4- e , 



(iH-e^)-— - -f(l-fei) 



80 



8ir 



(l+^:)^-(l+e3)^ 



(6) 



which we shall write as 1, m\, n\ ; and its final extension, to terms of the first order 



in t^', v\ to\ is 






1 • n 

ox 



In the same way we find that the other directions of principal strain in the final 
configuration are given by the direction-cosines 



and 



(14-61) ::^ 4-(l+62) 



dx 



dy 



(l + e,)^+(l4e3)f^^' 



dx 



dz 



{l + e,Y-{l+e,Y 



en/ . /. . \ dw^ 



■••? 



(l+^^3)-g-+(l+^2) g 

(14- 62)^- (1+63)^ 






dy 



cz 



{l-\-e,Y-{l-{-e,y 



1 



r 



(8) 



y 



(which we may write as 
in these directions are 



and 



■m\^ 1, n^2!> ^^^<^ -~^\j — n^25 l)? ai^d that the final extensions 






60 -1- 



3^' 1 






2 ' ^ ? 






63 + 



dz * 



(9) 



J 



The stresses in these directions, which we shall call the directions of x\ y\ and z\ 
referred to the original areas of the faces on which they act,^ are therefore 



Y' — Y -4- 9X^y '^^ , 9X^' 9^' , 9X^ 3i/ 



36i dx de. 



3eo ' S^J 



^3 



r^ ,»-t/JI» I 



20 

971 — 2 



/ . \ 3t^ . dv dw 
(m- 1 ) -TT- 4 -r- 4- TT- 

ex dy cz 



• • a . Qk/Kj* 



(10) 



* Of. the assumption of p. 193. 
2 O 2 
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Referred to the neiv areas of the faces on which they act, they are 









{l+e',){l+e',) 






71, 



71, 



(l+e'i)(l+e'2) 



and to the required degree of approximation we may write 



v/ ^ __ 



Jk. 



a; 



(l+(?2)(l+e3) „ 

2C 
m— 2 



1 



1 'dv' 1 3ic' 



1 + Cg 9?/ 1 + 63 " 3^ 



(l+e2)(l+e,,) 



(m— 1) 






3j/ 3t^/ 

9a? 3t/ 3^; _ 



s s t tt C3k/v^# 



• • • 



(11) 



Then if cc, ^/^ denote the co-ordinates in the final configuration, referred to the 
original axes, of the point which was originally at {x, y, z), so that 

we may find the stress components in the third configuration, referred to the original 

axes^ and to the strained areas of the faces upon ivhieh they act^ by the scheme of 
transformation 





x' 


y' 


z' 


X 


1 


—in\ 


-n\ 


y 


in\ 


1 


-t^2 


z 


n\ 


n'2 


1 



The following expressions are thus obtained (to the required order of approxima- 
tion) :~ __ _ -__ _ 



x_- 



m\{JJ^'-Y'y), 



m\ 






. » « • • l.L^I 



Now the stress-components (12) must satisfy the ordinary equations of equilibrium, 
which are three of the type 

• » » I X O / 



C^'V"- ^"Y"- P'Y"- 

— = — j =f — J- — =~ = U 

3a? 3t/ 3;^ 



• •«»«* 
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and since the co-ordinates of the point which ultimately goes to {x-\-Sx, y, z) were 
originally 

dX CX dx 



dx 

we have 

8 _ 1 8 dx 8_ dx 8 

9^"~ la.. j_§i^^^ (l+^i) {l+e,)dy {l+e^){l^e,)d~z' "*' 

1- "r c>2 "T _ 

OX 

It follows that (13) may be written (to our approximation) as follows : — 



{l+e,){l+e,){l + es) 



1 + ^2 3a; By 1 + e., 3« 3a3. 



0. . . (14) 



+ — < (m — 1 ) — - -\ j 

m — 2 I dx^ dxdy dzdx. 

, {l+e^)X^-{l-\-e2)Y ^dm\ ^ {l-{-e^)X^-{l-{-e^)Z,dn^^ 

~\~ ' — ~"^ —' — ' — ~i~ ' 

1+62 8^ 1+^.3 dz . 

Substituting for m\, n\, we have finally 

m— 2 8a^^ 8^^ 8^;^ m— 2\8,x8t/ 82:83^/ 

+ X.+Y, / 8V ____ _8V_\ ^ X.+Z, / 3y _ ^ ^ (15) 

iCil 1 ^i + ^2 \ \9^' dxdy) ^q/i , ^i + ^3 \ \82;' 82^80^/ ^- • V / 

X 

and two similar equations. In any ordinary problem we may neglect ei-~ ... in 

X 

comparison with 7=f . . . . 

The equations thus obtained may also be written (with Lame's notation for the 
elastic constants) as follows :• — 

where 

. / __ 8u' dv^ dw^ 
dx dy dz 

dx^ dy^ dz^ 

2 ' = ^' — -^ 2 ' = — — -^ o ' — ^ _ ^ 
"" 8^ 82; ^ 82; 3a? "^ Sir 3^ ' 
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and in this form they may be conveniently compared with the ordinary equations of 
elasticity. "^ 

The three equations of the type (l 5) we shall term Equations of Neutral Equilibrium, 
The equilibrium of the stress-system X,,, Y^, 7a^ will be neutral, provided that solutions 
for %i\ v\ %v^ exist which satisfy certain boundary conditions. These boundary 
conditions are peculiar to each problem, but usually express the condition that the 
additional stresses involved by u\ v\ iv^ shall vanish on certain boundary surfaces. 
They never deterwiine the magnitude of u\ v\ w\ so that our solution gives the form 
only of the distortion which tends to occur in the body under consideration when its 
equilibrium becomes unstable. It gives a definite relation between the stress-system 
X^... and the dimensions of the body, which must be satisfied in order that any 
distortion may be permanent ; but if this relation be satisfied, no limits are imposed 
by the equations upon the magnitude of the distortion which jnaj occur, f 

Example in Rectangular Co-ordinates, Stability of Thin Plating under 

Edge Thrust. 

It seems advisable, before we employ a new method on problems which have not as 
yet received satisfactory treatment, in some degree to test its validity by the result 
to which it leads in a more familiar example. For this purpose we may consider the 
stability of an infinite strip of flat plating under edge thrusts in its plane. The 
accepted formula^ for the thrust necessary to produce instability, per unit length of 
edge, is 

m^—i r 
where 

2t = thickness of plate, 



j5 — 3 2 — 7 '^^ 7i ' ' ' * ' * \ ■'• ' 7 



/ = breadth of plate, 

and the opposite edges are simply supported. If the edges are built in, the thrust 
required has four times this value. 

To investigate this problem by the new method we take axes Ox and Oz in the 
middle surface of the plate, in the direction of its breadth and length respectively, and 
Oy perpendicular to the middle surface. The initial stress-system is then given by 



X^ = const. = G (say), 
Y — 7 = • 



> ■ / 1 R ^ 



^ Love, op, cit., § 91, equation (19). 

t The equations are, however, rigorous only in the case of infinitesimal displacements; cf, footnote, 
p. 240. 

X Cf, Love, op, cit., § 337 (a-), whence the above expression may be obtained. 
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and we may assume that the system of strain which is introduced at collapse will be 
two-dimensional, so that 



dz dz 



0, 



z 



— const. 



• • • • • « • \X %7 J 



The third equation of neutral stability (for the direction Oz) is then satisfied 
identically, and the other two equations become ( if we neglect terms of order j^^u^ -** 



\2^J 



ni 



?h/ 



G a ?a/ di/ 



ano. 



2 dx^ ' dy^ on — 2 dxdy -kGdyKdx 



ni' 



3V o ^^— 1 8^^?/ 

dx 



m 



dhi^ 



'JjOlJ 



%. 



G 3 /du^ dl/' 



m — 2 dy^ m-~2 dxdy ACdx\dy dx, 



^ 



0, 



V 



• • • \ -l^ \J f 



0. 



J 



Let us assume a solution of the form 



u 



2 [Ua sin a {x-\-X)^)^ 



i/ ^ S [V« COS a {x + .Xo)] 



(21) 



where IJ« and Y^, are functions of y onl)^. It is easy to show that this assumption as 
to the phase-relation of u^ and i/ is justified. We have then 



2 — — — -a'^U«+ 1+ ^ 



m 



m — 2 



G\ c?V 






40/ (^?/ W~2 40/ c?^ 
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m 



40. 
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v^Tl 






— 2 40/ c?^ 



0, 



"\ 



V 



0. 



(22) 



The solution of these equations is of the form 

U^ = (P|/ + Q) sinh a^ + (R|/ + S) cosh ay^ 



r' 



V,= --< R^+S 



' 3m-~4 G 
m — 2 



V-. 



m 



-2 40 




> sinh «?/ 



\ Ps/+Q 



iim—4: G \ 
m-2 4C\R 



1 



/- 



I 



m 



'???. G la 

2'~4C/ 



>- 



cosh ay, 



J 



-J 



(23) 



where P, Q, R, S are constants. 

The boundary conditions now demand attention. It is clear that the stresses 
introduced by %i\ v\ tv^ must vanish at the surfaces of the plate. Hence these 
surfaces will still be planes of principal stress, and, moreover, the normal stress upon 
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them must vanish. But, as we have akeady seen, the Une which becomes a direction 
Oy^ of principal stress has initially the direction-cosines 

-mi, 1, n^: 

it follows that at the surfaces of the plate the expression for mi and % must vanish 
identically ; moreover, at these surfaces, Y^^' must vanish. These conditions may be 
written in the form 



dy 



(1+62) — + (1 + 63) 



^ 



0, 



(l+eO^ + (l+e,)^=0, 



dy 



dx 



> identically, when y = ±t. . . (24) 



(,n-l)— +-^+_=:0, 



J 



The first condition is already satisfied. The other two give ( if we neglect terms of 



order —u^ ... 



and 



(1+eO ^-(l+e,)«V„ = 0, 
(m-l)^+«U, = 0, 



^ 



> when y = ±t^ 



dy 



» » 6 



. . (25) 
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and 
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G^ 
40 



m 



w 



+ R 



hen 



G 

m^^iC 

G 



cosh ay 



+ Q (m— 2) a sinh ay 



(m~2) ay cosh ay— 2 (m—l) 



1 + 



m. 



40 



_G 

m-2 40 



y = it 



sinh ay 
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(27) 
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Thus we obtain 
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There are two solutions of the equations (28). Either 
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— at cotha^, 



. . \Zl J J 



at tanh at. 



and 



Q = R -= 0. 



• s • 



e • • 



. . . (30) 



The criterion for neutral stability is in the first case 

2m^— 1 



at (coth otf. — tanh at) 



and in the second case 



at (tanh. a^ — coth at) = 



1 + 



1 + 
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so that the values of G, for which collapse by mstability may be expected to occur, 

are given by 

G __ l~-2a^ cosech 2at /^^x 

2m(m + l)C 2m^—l — 2a^ cosech 2a^ 
and 

G _ l+2a^ cosech 2at 

277^(m+l)C 2m^~l + 2a^ cosech 2a^ 

respectively, the total thrust, per unit length of edge, being 



. • . . . . \oZij 



Jl? — ' ' /it\JJ (dt)y 

The first approxinmtions to a solution, in terms of t, are 



and 



G- 


3 • 1 

m— 1 


X — 


gm+l^ 



: -E ........ . (35) 



respectively. Since the complete wave-length of the corrugations into which the 
plate distorts is 

f\ — ,, « « • . . • . , . , \ tj\j t 

we see that (34) is equivalent to (17), and that the latter formula is therefore 
supported by our investigation as a first approximation. The second solution (35) is 
without practical interest, owing to the magnitude of the thrust required to produce 
collapse. It refers to a type of distortion, theoretically possible for an ideal material 
without limits of elasticity, which is approximately realized in actual specimens of 
ductile material, when tested to failure under compressive stress. Since Q = R = 0, 
we see from (23) that in this type the middle surface remains plane. In the first 
type of failure, where P = S == 0, we find that U, ^ when ?/ = 0, so that the middle 
surface undergoes no change of extension in the distortion given by te', v\ tv'f^ 

Equations of Neutral Equilibrium in Cylindrical Co-ordinates. 

Deprivation of the Equations, 

The equations (15) of neutral equilibrium are expressed in a form which is 
unsuitable for the investigation of problems concerned with the stability of thin 
tubes, and we have next to obtain the corresponding equations in cylindrical 

* Besides the harmonic sohitions to (20) we may have 
but g, h, and k vanish in virtue of the boundary conditions. 
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co-ordinates. We shall limit our discussion to stress-systems which produce a 
displacement symmetrical about an axis, up to the instant at which the equilibrium 
becomes unstable and distortion occurs : in Pearson^s notation, the principal stresses 

in the equilibrium configuration are m 00^ and zz, and these quantities are functions 
of r only. 

The new equations are derived by a method very similar to that which has already 
been explained. The co-ordinates of a point in the unstrained configuration are 

n &, ^; 

in the second configuration (of equilibrium) they are 



and in the third configuration (of slight distortion from the position of equilibrium) 
they are 



r-hii-hu^ 



+ 



V 



r+ti 



.z + tv+w^ 



(the radial, tangential, and axial displacements ii\ v\ tv^ being ultimately taken as 
infinitesimal). 

The extension of a line-element joining the point (r, 0, z) to the point (r+^r, d-}-S6, 
z + Sz) is 



(/ = —1-1- 



a/I 



l-^-fTi^-j^n^Jl 



dr dr \r dd rj dz _ 






dr 



4-m (1-1 1 ^-~T7r) +^^ "T~ 

(jZ 



r T r ad 



where 












r de 



dz 



dz 



• • • • 



m = r 



S9 



and 



71 






(37) 



and this has a stationary value for a line very slightly inclined to the radius, given 
by 



(l+^i) 



m, = 



\r 00 r I or 

cm . /-, . \ (hv 



(l+e,)?^ + (l+e,)^ 



m 



dz 



dr 



{1 + e,f - {I + e,Y 



. • . . . \0<J ) 



J 



where 6^, eg, e-s are written for dti/cr, u/r^ and dtv/dz respectively, 

2 D 2 
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The extension of this hne-element in the third configuration is 



e/ = ey 4- 



cm 



e • fr • e 



. . . (39) 



and its inclination to lines issuing from the point (in its final position) in the radial, 
tangential, and axial directions is given by the direction cosines 



1, m\, n\, 



wnere 



(1+6,) 



'1 37/ 



r]%\ 



T dO 



V 

r 
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OZ ■ f. 



dr 



(1 +63)^- (1+63)' 



• 9 
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(40) 



We find also that the other directions of principal stress in the final configuration 
are initially inclined to radial, tangential, and axial lines through (r, 0, z) at angles 
whose direction cosines are 



a 



jid 



m,, 1, n.,, 



'^^•'15 "^"^^2? '- 



and that in the third configuration they have corresponding inclinations to radial, 
tangential, and axial lines through (r, 0, 2;)-~"in its final position-™~which are given by 



771;^, 1, n 2? 



jid 
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n ^ ~ 
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The extensions of the corresponding line-elements in the final configuration are 
respectively 



1 at/ 



/ u 

r r a 



2 



and 



e'. 



e-A -4- 



00 



» • 



(42) 
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so that the principal stresses in the third configuration, referred to the final areas of 
the faces on which they act, are 



r^nr 



/y^/y 



Tl 
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Then by the scheme of transformation 
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• (43) 



we find for the stress components in the third configuration, at the point which was 
initially at (r, 0, z), — referred to axes in the radial, tangential, and axial directions 
through the final position of the point, and to the final areas of the faces upon which 
they act — the expressions 



^ = 07^ = m\ {rr'-ee') 



ee = ee', zz = zz', 



m\ 



[( 1 + eO rr~{ 1 + e.,) 66] (to the 



(l+ei)(l+e2)(l+e3) 
required approximation), &c. (44) 
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The equations of equilibrium, to be satisfied by the stress components (44), are* 
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the equations (45) may be expressed in differentials with respect to r, 0, z. The 
terms which do not involve ii\ v\ w' vanish in virtue of the equilibrium conditions 
for the second configuration, and only tei^ms of the first order in these quantities 
need be retained. 

In general, e^, e^, <h i^^y ^ be functions of r, but in this paper we shall only 

consider problems in which e^ and zz have constant values. Moreover, in all problems 



2 



of practical importance we may neglect terms of order -7^?/... in comparison with 



■^ Love, oj), at,, § 59 (i). 
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terms of order jq- id — We then obtain, as the eqnations of neutral equilibrium in 
cylindrical co-ordinates, 
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(48) 



Equations (46-48) represent the conditions for neutral stability in the equilibrium 

of a body subjected to a stress-system rr, 66, zz, where zz is constant, and rr and 66 
are functions of r only, which satisfy the condition of equilibrium 



orr rr—66 



cr 
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r 



(49) 



For comparison with the ordinary equations ot equilibrium in cylindrical 
co-ordinates they may be written (with Lame's notation for the elastic constants) in 
the forms'^' 
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* Cf, Love, t^jx ciL, § 199. 
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Examples in Cylindrical Co-ordinates, Stahility of Boiler Flues and 

Tnhular Struts. 

The Equations of Neutral Equilibrium in Cylindrical Co-ordinates enable us to deal 
successfully with some difficult problems connected with the stability of cylindrical 
tubes. Two examples of considerable importance will be discussed in this paper— the 
collapse of boiler flues and the strength of tubular struts. It should be noticed that 
neither of these problems has been quite satisfactorily treated by the ordinary theory 









Fig. 2. 



of thin shells, which requires the assumptions that the middle surface of the shell is 
unextended, and the inner and outer surfaces free from applied tractions'^ ; hence their 
solution is a problem of considerable interest, even apart from practical considerations, 
and has attracted a great deal of attention. It will be convenient at this point to 
review the work which has already been done. 

The question of the stability of tubular struts is important, owing to the frequency 
of their employment in practice. In economy of material the cylindrical tube 
possesses an advantage over struts of solid cross-section, and both the theory of 



^ Cf. A. B. Basset, ' Phil. Trans. Roy. Soc./ A, CLXXXL, p. 437 ,; and Eayleigh, ' London Math, 
Soc. Proc.; vol. XX., p. 379. 
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EuLER^ and LAGRANGEt and the more practical formula of Eankine suggest that 
this advantage increases without limit as the thickness of the tube is reduced. Such 
a conclusion is, however, inaccurate, for types of distortion are possible in the case of 
a tube which do not involve flexure of the axis, and when the tube is thin these 
types, of which some practical examples are shown in fig. 2, may be maintained by a 
smaller thrust than would be required to produce failure of the kind discussed by 
EuLER. Moreover, the natural wave-length, for these symmetrical types of distortion, 
is in general small, so that distortion can occur without hindrance in quite short 
tubes. Hence, for a considerable range of length the strength of a tube to resist 
end-thrust is practically constant, and is not given by any of the usual formulae for 
struts. 

The determination of the strength of tubes to resist these symmetrical types of 
distortion is obviously a problem of the highest practical importance, and has 
attracted a great deal of attention in recent years. Illustrations of collapsed tubes, 
showing symmetrical types of distortion, have been published by A. MallockJ and 
R. LoRENZ,§ and a great deal of experimental work has been carried out by 
W. E. Lilly. II Theoretical discussions, by approximate methods, have been proposed 
by A. GROS,f W. E. Lilly,'^^ S. TiMOSCHENKott and R. Lorenz.JJ 

The problem of the boiler flue seems first to have been suggested by the experi- 
ments commenced by Fairbairn in 1858. §§ These showed that the collapse of tubes 
under external pressure was in some degree analogous to that of straight columns 
under end-thrust, and a discussion of the phenomenon, based on Euler's theory of 
struts, was given by W. C Unwin,!! who assisted Fairbairn in his research. The 
similar problem of a circular wire ring subjected to radial pressure has been 
discussed by M. BresseIFIF and M. Levy,'^^'^ and rational theories of the boiler- 
flue problem have been given by G. H. BRYAN,ttt A. Foppl,JJJ P. Forghheimer§§§ 

^ " Sur la force des colonnes," * Hist. Acad. Berlin,' XIIL (1757), p. 252. 
t " Sur la figure des colonnes," ' Miscellanea Taurinensia,' V. (1773). 
t ' Eoy. Soc. Proc.,' A, LXXXL (1908), p. 389. 
§ * Physikalische Zeitschrift,' XII. (1911), p. 241. 

II ' Proc. Inst. Mech. Eng.,' 1905 ; 'Inst. Civ. Eng. Ireland,' 1906 ; ' Engineering,' January 10, 1908. 
51 ' Comptes Rendus,' CXXXIV. (1902), p. 1041. 
^^ 'Inst. Civ. Eng. Ireland,' 1906. 

ft 'Zeitschrift f. Mathematik u. Physik,' LVIII. (1910), p. 337. 

XI ' Zeitschrift des Vereines Deutscher Ingenieure,' October 24, 1908 ; 'Physikalische Zeitschrift,' XII. 
(1911), p. 241. 
§§ ' Phil. Trans. Roy. Soc.,' CXLVIIL, p. 389. 
nil 'Proc. Inst. C.E.,' XLYI. (1875), p. 225. 
^H ' Cours de Mecanique Appliquee,' 1. Partie, Paris, 1859. 
**^ ' Liouville's Journal,' X. (1884), p. 5. 
ttt 'Proc. Camb. Phil. Soc.,' VI. (1888), p. 287. 
tJt 'Resistance des Mat^riaux ' (1901), p. 286. 
§§§ ' Zeitschrift des Oesterreichischen Ingenieur- und Architekten-Yereines,' 1904. 
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and R. Lorenz.^ W. E. LiLLYt has indicated the correct form of the result for 
an infinitely long flue, and A. E. H. LoyeJ has discussed the strengthening effects of 
constraints which keep the tube circular at its ends. 

A. B. Basset§ has given a very clear exposition of the difficulties which are 
encountered in an attempt to construct a theory of flue collapse by usual methods. 
To obtain sufficient equations we must assume that the middle surface undergoes no 
extension ; and the existence of pressure on one or both surfaces of the tube not only 
makes this assumption very improbable, but violates an essential condition upon 
which the theory of thin shells is based. When one surface only is subjected to 
pressure, there is reason to believe that Bryan's solution is substantially correct ; 
but no treatment can be looked upon as rigorous which neglects the cross-stresses in 
the material. 

The experimental researches of A. P. Carman|| and R. T. Stewart^ have revived 
interest in this problem, since they offer the first information which has been obtained 
as to the behaviour under practical conditions of tubes which in circularity, uniformity 
of thickness and homogeneity are fair approximations to the ideal tube of theoretical 
analysis. '^'^ 

We commence our discussion by considering the stability of a thin cylindrical tube, 
subjected to the combined action of end and surface pressures. We shall thus be 
able to derive the required solutions for the thin tubular strut, and for a boiler flue 
without end-thrust, as particular cases, and from the general solution we may obtain 
indications of the way in which end-thrust tends to promote the collapse of a 
boiler flue. 

In the most general form of the boiler-flue problem, as enunciated by Basset, ft 
pressures are acting on both surfaces of the tube, and we shall therefore investigate 
conditions for neutral stability in a tube subjected to the following system of 
stresses :~ — 

(i.) An end- thrust of total amount ,S, uniformly distributed ; 
(ii.) An external hydrostatic pressure, of intensity ^j ; and 
(iii.) An internal hydrostatic pressure, of intensity 1^2- 

^ ^ Physikalische Zeitsclirift,' XII. (1911), p. 241. 

t 'Inst. Civ. Eng. Ireland/ 1910. 

} ' Proc. Lond. Math. Soc.,' XXIY. (1893), p. 208. 

§ ' Phil. Mag.,' XXXIV. (1892), p. 221. 

II '^Eesistance of Tubes to Collapse," 'Bulletin of the Univ. of Illinois/ No. 17, 1906. 

^ " Collapsing Pressures of Bessemer Steel Lap- Welded Tubes," ' Trans. American Soc. Mech. Eng./ 
1906, p. 730. 

^"^ The experiments of Fairbairn were restricted to tubes which were constructed from sheet metal, 
with brazed and riveted seams. 

tt Log. city p. 223. 



ME. E. V. SOUTHWELL ON THE GENEEAL THEOEY OF ELASTIC STABILITY. 211 



We shall consider a tube of indefinite length, of which the inner and outer radii are 



and we shall write 



a±t (so that the thickness is 2t)^ 



T for the ratio -• 

a 



The corresponding stress-system, for the position of equilibrium, is easily obtained."^ 
We have 



n^/T/» 7—;; __ 



1 



a' 



^.(l+r)^-f.(l-rN^(ii-ft)(l-T^) 



ee = - 



1 



r 



a' 



^ 



^x(l+T)^-f3(l-T)^+3(fl-f2)(l-T^) 



r 



and 



zz 



% 



4:irat 



. . . (51) 



J 



It can also be shown that 63 is constant, and equations (46-48) may therefore be 

taken to express the conditions of neutral equilibrium. The degree of approximation 

to which these equations have been obtained (p. 206) will be maintained for the rest 

2 
of this paper, i.e., terms of order -^ ^6^.. will be neglected. They may also be written 

as follows : — 



^m— 1 /3V 1 du^ 



U 



f\ 



m-— 2 \ 3r^ r 9r 7' 



i, + (i + AU9!!f: 



+ \1 + 



A 

4 



1 + 



r^ 4 J dz 



a 



+ 



+ 



_ A\ 1 9"^' _ / 3m -4 , A\ 1 ^ 
.m-2 2/^3^30 \?n-2 2 / r^ 80 



m 



A 



m—2 



1 + 



(ra"\ . Bl dhv' 



+ 



4 J dzdr 



= 0, 



(52) 



m, 



^m 



_ A\ 1 9V , / 3m -4 A\ 1 ^ 
-2 2/ r 3^30 \?n-2 '2/r'30 



+ 1 + 



"T" 1 -L I 



A\ /av . 1 3^^' ^' 



2 
A 






o m— 1 j^ 3V 
^^ m — 2 r^ dff^ 



1 



n^2 



Bl 3V f 7n 
+ 



3^^ 



m—2 



A /' __ era 



r 



'\ . Bl 1_3V 
4 J r 30 3^; 



-f 



0, (53) 



and 



m 



A 



\m—2 4 



1 + 






'\ . Bl 3V 



+ 



+ 



m 



A 



+ 



m 



m — 2 
A 



^^1-^^) + 



4 J 82; 3r ' l,m— 2 
^^ Bl 1 3V 



1 






^\ . Bl 1 du' 



+ 



r 



2 



4 J r 30 3^; 



, ,, , A / , , o-a 



2^ 



ry%'^ 



r dz 

Bl 3V 

3r^ 



+ U + 



1- 



o-a' 



2^ 



ry%' 



,2 



4 1 r 3r 1 4 \ r 



2 



Bl 1 3W , o^— 1 9^^' A /^^\ 



4 r^ 30 



m — 2 dz" 



* Love, ' Mathematical Theory of Elasticity \ (2nd edition), § 100. 

2 E 2 
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where 



A= -^[%{l+rr-%{l-rn 



B = 



^ 



AiratG 



'^ t • t e > * • \ KJ KJ J 



and 



cr 



^l(l+T)^-i2(l-T)^' 



J 



We may assume a solution for equations (52-54) of the form^ 



%i 



V 



IV 









2 



U;,, q sin h{Q-^ Qj) sin ^ (^ + 0^) 

V;, „ COS /c (0 + 0;;.) sin 2 (^ + ^ ) 

W;,, ^ sin /i: (0 + 0^.) COS " {z + z^ 



^ 



> 



ft S 9 



^ 



• (56) 



where ^ must be integral, and U;,^^, V;,,^, W;,^^ are functions of r only, which satisfy 
the differential equations 



9 ^'^1 ( ^^ Id 1 
m— 2 Vdr^ r dr r^ 

m 



F/. . A \ gM. . A /, . o-a^ 



r^ \ 2 



-^ 1 + 



a 



k 

1 

a 



jn—2 2 J r d7 



AM c? /3m-4 , A\ 1 

+ 



m — 2 



A /-, (ra^\ B 
I H — -- ) + — 

r^ / 4 



m~2 ' 2 / r^ 
d 



4 

V 



1 + 



r 



k,q 



dr 



^^ k,q U, . 



^1 
4 J 



u., 



• » • 



(57) 



h 



i} 



on 



,m— 2 2 r dr \m—2 2 J r^ 



u 



^>g 



4 



r 



i + AV-i + i^ 



2 / \dr^ r dr 7 



Y-^j on— -2 r^ a? \ 



A 



1 



era' 



k1 

a 



and 



m 



A 



m — 2 4 



1 



o-a" \ B 



r 



4 






0. 



• • • 



a 



{ on 



-—(1 + 

m~2 4 \ r\. 



(ra"\ . Bl d 



+ 



k^ 

a 



on 



4 J dr 

,2^ 



+ 



r 



m 



A /. 



o-a' 



m — 2 



-f 



B 



1 

r 



U 



B 

4 



V 



/eg 



• • e • 



^. g 



m — 2 



A /-J __ q-a^ \ B 



r V ■ 4 



ly 



A, g- 



+ 



i + A(n_-^ 



4 



/.2 



(ir^ I 4 \ 



A /. _^ era 



2^ 



/yt^ 



,2 



T dr 



r 



. (58) 



2 ^^lui . 2!" 
m—2 a^ 



W,,, = 0. , 



• • 



(59) 



It is easy to show that the phase-relations assumed in equations (56) are necessary. 
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The boundary conditions now require investigation. From the consideration that 
the cyhndrical boundary surfaces of the tube must continue to be tangent to principal 
planes of stress, in any possible type of distortion, we deduce the conditions 



m. 







, identically, when r = a±t. 



(60) 



Ui ~ 



The other boundary conditions are more complex. Since the pressures acting on 
the surfaces of the tube are hydrostatic, it is clear that the radial stress, as defined 
on p, 193, is increased at points on the boundary surfaces of the tube tvhere the 
distortion involves positive extension. In the notation employed above, we have 



rr^ = —Pi, when r = a + t^ 
= —^2? when r = a—t, 
and from (43) we deduce the following equations, which must be satisfied identically, 



# 



2 



m— 2 



(m-1) 3™ + " + ~ -^ + -^- 

dr r r ad cz 



•\ 



r 



\. 



ii 

c 
c 



u' . \^v' dw' 



r r do 



+ 



dz 



?/ 1 3;?/ dw^ 
r r do dz 



when r =^ a-\-t, > 



5? 



r = a—t, 



J 



(61) 



Substituting from (56) in the identities (60) and (61), we finally obtain, as the 
required boundary conditions in U;^,^, ^k,qi ^^d W^^^^, 



^U,„+{(l-A.^)|-^^V,„ 



r 



1 

r 



^ 



= 0, 



2U..,+ 1 



A / , , a 



Bl d 



•V 2 dr "■" 



W. , = 0, 



t 



and 



2^^=^iu,,-f 



2 



m 



2 dr 



m—2 



A 1 + ^ 



aa^] n jj k y ^_ 2 w 






T^ ; 1 



k,q ^ k,q ^ ^ k,q 

r r a 



0. 



^ 



(62) 



when r = a-\-t. 



^ In obtaining these equations it should be noticed that before distortion occurs - JPi and - 3^2 are the 
values at the boundary of 



n-'/y* 



(l+^2)(l+^3)' 



and not of rr, if we retain the significance for rr which was assumed on p. 193. The distinction is not 
really needed for the approximation of the following work, but it may lead to confusion if neglected. 
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The differential relations (57-59), with the boundary conditions (62), are 
theoretically suflficient for an exact solution of our problem : we shall, however, 
content ourselves with approximate solutions for (Hi — IP2) and %, correct to terms 
in T^. To obtain these, we assume solutions for U;,, ^, V^^ ^, W;,, ^, in series of ascending 
powers of the quantity {r—a). Thus we write 



<? I e 



'^k,q = fo + fl 






S 



2 A? 



\ 



1 s^ 



2! 



.a 



V 



hq 






\^^/ 



\a/ 2l\a/ 



» • 



. . . (63) 



y 



where r = a+h. 

We may now derive, from equations (57-59), any required number of relations 
between the undetermined coefficients ^o="?7o---fo--^ and the boundary conditions (62) 
take the form of equations in series of ascending powers of the small quantity t, in 
which the sums of the odd and of the even powers must vanish separately. If we 
neglect in these equations terms of order higher than some definite power of r, we 
may obtain corresponding approximations to the values of A and B, by the elimination 
of the undetermined coefficients. 

The approximate boundary conditions, correct to terms in r^ are^' 



/cfo + ^Y^2-';o+{l-A(a— T')};7i + (i-A)T';;2 + (l+A)y^3 = 0, 



(64) 



^a + ^^-j^3+A(o---"T^)^i+{l~~A(o---T^)}^,+ U— . . (65) 



2 



qio+qY^^+ 



^ - ■ I 1 _ A(l+(r) + B _^ ^^ ^^ _ ^^,^^ ^_ f J B \ T 



2 



2 , 2 ^3 - U, . . . 



• ((56) 



Q'fi + gy6 + ^('^-2T')^i + jl- 



[. A(l+q-) + B , 3 



B\ a 



2 



+ iAr^^r f.-iArYa + i 1-^)^^4 = 0, (67) 



1 _™_2 A(l +0-3x^)1 |o+ {m-l-(m-2) Ar^} ^, + (2m--l) ^ ^3+(ot-1) ^ i, 



q\l-'^A{l + ^-r')\^,-q{l-{m~-2)A}T%-q^^, 



0, . , . . . (68) 



* In deriYing these boundary conditions it is to be noticed that o- is to a first approximation aqua 
to - 1, so that to our approximation - t^ may be written for o-t^. 
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and 

(to-2) a {<T-2r') fo+ [m-"^ A (1 + ^-3t^)| a+ jm-1- ^^ At^}- ^, 



_2 



+ (3m-2) ^ i^+{m-l) ^ i^-{m-2) IcA {a-2r'),,,-Jc\ I-^^ A(l+o— 3t^) 



6 



6 



2 



'?! 



_2 



+ (m-2)M^,3-/«y%-g|l-^A(l-<.+r^)Ko 



g^l_^A(l+a-T^)|A-g{l-(m-2)A}^f,-g^f3 = 0, . . (69) 



2 



2 



6 



and these equations involve the fifteen coefficients 

fil • • • t,4^ % • • • '?4) So • • • M- 

From equations (57-59) we may obtain nine other relations between these 
coefficients, as follows : — 



2 
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(78) 



We may now eliminate the coefficients from equations (64-78), and obtain a 
determinantal equation, of fifteen rows, which gives a relation between A, B and the 
dimensions of the tube. This relation is the condition for neutral equilibrium of the 
initial stress-system, and is clearly correct to terms in t^ ; but by further consideration 
of the terms inyolved we may show that the labour which w^ould be required for its 
complete evaluation is unnecessary, and as the fifteen-row determinant may be 
written down directly from the above equations it will not be given here. 



Sohction for Boiler Flue without End Thrust 

We shall begin by deriving a sufficiently approximate expression for the difference 
of pressure required to produce collapse of a tube, when there is no resultant end 
thrust or tension ; and in the first case we shall deal with a form of collapse possible 
only in the case of tubes of infinite length. That is to say, we make B and q equal 
to zero in the fifteen-row determinant, which may then be reduced to one of ten 
rows. 

In the latter determinant we may treat A as a quantity of order t^ ; for if A be 
put equal to zero, and the determinant be expanded, the terms which are independent 
of T vanish identically. Hence —A may be written for a-A, and Ar^ may be 
neglected. 

The ten-row determinant, simplified by these and other obvious modifications, is 
given on pp, 218 and 219. Expanding it from the top row, with the neglect of terms 
of order higher than t^, we obtain 
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But to a first approximation 



A=-^(fx-f,), 



SO that we have 



m 



m— 1 



,2 //- \r> 



'A 



'"^^ E(F-1)(^~), . „ (80) 



which agrees with Bryan's result/^ 

To complete our discussion of this problem we must consider types of distortion in 
which the axial wave-length is finite, and thus obtain a theoretical estimate of the 
strength of short flues with fixed ends. A solution giving A correctly to terms in t^ 
may be derived from the complete fifteen-row determinant ; but we may show that 
for practical purposes the labour which this evaluation would entail is quite 
unnecessary. 

We find first of all that those terms in the expression for A which are independent 
of T contain q^ as a factor. Now 2AC being approximately equal to the mean hoop 
stress in the tube before collapse, it is clear that A must in all cases of practical 
importance be a very small quantity. It follows that in the expanded equation the 
terms in A are of primary importance, and A^ and higher powers may be neglected ; 
further, since q must also be small, that terms in q^ and higher powers of q may 
be neglected in comparison with terms in g'^, and that of the terms in r^ those which 
involve q are negligible in comparison with the terms already found. 

In accordance with these principles we may derive the terms which are required to 
complete our solution from a nine-row determinant, obtained by omitting terms in r^ 
from the general determinant. This simplified determinant is given on pp. 218 
and 219. Further, we may neglect A^ in the expansion, and in the coefficient of A 
retain only those terms which do not involve q ; we thus obtain the equation 



But, by equations (55), 






4Ct 



and therefore, to the approximation of equation (81), 



fx-f.-4^C^|-^.--2E-j^---^^ ^(82) 



-H- 



Cf. footnote, p. 209. 
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Combining this result with (80) we have, as our final expression for the pressure- 
difference which can produce collapse of the flue, 



%-}% = 2E ^ 

a 



<t 



+ 



mf 



k'{k'-l) ^m?~l 



(F-1) 



OJ 



2 



(83) 



In this equation tja is the ratio of the thickness to the diameter of the tube, and h 
is the number of lobes in the distorted form of its cross-section. The quantity q is 
connected with the axial wave-length X of the distortion by the relation 



qX = 27ra. 



(84) 



We may imagine a flue subjected at its ends to constraints which merely keep the 
ends circular, without imposing any other restrictions upon the type of distortion."^' 
In this case the end conditions may be written in the form 



It' - 
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> when z — < 
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(85) 



and from (56) it is clear that Z, the length of the flue, is equal to irajq. 
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P Added June 8, — Thin circular discs, inserted into the tube at its ends, but not fixed to it, would 
approximately realize these conditions.] 
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In practice, the end constraints will also tend to maintain the cylindrical form at 
the ends of the flue, and this effect will strengthen the tube, by an amount which is 
not easy to determine exactly. In any case we may say that 

I 1 

a q 

and we may illustrate the way in which the end effects die out by plotting the 
pressure diflPerences QPi— P2) against the quantity g""\ To do this we must take 
some definite value of the ratio tja^ and plot different curves for the values 2, 3, ..., 
&c., of h The result is shown by fig. 3, in which the following values have been 
assumed for the constants : — ■ 

E = 3x10^ pounds per sq. inch, 
:^ 2*07 X 10^^ dynes per sq. cm. 

T Q If "I 

a 

From an inspection of the different curves we see that long tubes will always tend 
to collapse into the two-lobed form, since the curve for h -=2 then gives the least value 
for the collapsing pressure, but that at a length corresponding to the point A the 
three-lobed distortion becomes natural to the tube, and for shorter lengths still, of 
which the point B gives the upper limit, the four-lobed form requires least pressure 
for its maintenance. Thus the true curve connecting pressure and length is the 
discontinuous curve CBAE, shown in the diagram by a thickened line. 

Whatever be the relation between q and the length of the flue, it is clear that 
instability is theoretically possible in cases where the distortion involved is not even 
approximately '' inextensional/' For if r is sufiiciently small, the collapsing pressure, 
as given by (83), need not involve elastic break-down in the position of equilibrium, 
even though the first (or '' extensional") term in (83) be equal to, or even greater 
than, the second. Of course, elastic break-down will occur by reason of the extension 
very soon after the commencement of the distortion. Nevertheless, failure in such a 
case must be i^egarded as due entirely to instability ; for if this source of weakness 
were removed, effective resistance could be offered for an indefinite period to pressures 
which actually result in collapse. 

Comparison with Experimental Results, 

Although, as we have just remarked, it is theoretically possible for failure to occur 
by true elastic instability in comparatively short tubes, yet the relative dimensions of 
the tubes must be such as it would be quite impossible to test experimentally. In 
any practical case, instability will not occur until the properties of the material have 
been altered by overstrain, and the value of the pressure at collapse is therefore very 
much less than the foregoing theory would suggest. 
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It is, however, of interest to compare the general shape of the theoretical curve 
CBAE (fig. 3) with the results of experiment, and fig. 4 has been constructed for this 
purpose. It represents a number of tests conducted by the author upon seamless 
steel tube (0*028 inches thick and 1 inch in external diameter), and shows the relative 
amounts of resistance to external pressure offered by diiferent lengths of tube. In 
these experiments (selected for fig. 4 from a more comprehensive series which is still 



1 



^8 



h 






5 

.S o 
^ o 



^ 



.5v 



o 
o 



^ 



^ 

^ 















•Oo 
o 
















Steel Tubes. 
(Tested by hydraulic pressure.) 
) lobes in distorted cross-section. 

^^ )3 L» >JJJ 
** » J5 » 







4> Twc 






t 


00 Tlir( 
s^-^ Foil 




t 














•i 


^ „_ - :.__.. 














^ 
^ 




■¥ 












- 













UneupporteiL 5 Lengtih of 10 Tube in 15 IivcJies. 20 

Fig. 4. 



in progress) the ends of the tube were gripped by means of slightly conical plugs and 
sockets, the interior being kept in free connection with the atmosphere, and no 
attempt was made to balance the axial thrust due to hydrostatic pressure on the 
plugged end of the tube. Other experiments have shown that the existence of this 
thrust is not seriously important. 

It will be seen that the general shape of the theoretical curve is well reproduced, 
as well as the changes in the number of lobes which characterize the distorted cross- 
section. Similar results to those of fig. 4 have been obtained by Carman,^ but his 
experiments were not sufiiciently numerous for a satisfactory comparison with the 
theoretical curve of fig. 3, his object in conducting them being merely to discover 
what is the limit of length beyond which the strength of a tube may be taken as 



■5f 



Cf, footnote, p. 210. 
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sensibly the same for all lengths. The main interest both of Carman and of 
Stewart"^ was confined to tubes in excess of this limit, experiments on which may 
fairly be compared with the theoretical formula (80) ; their results showed that this 
formula gives a satisfactory estimate of the strength of very thin brass and steel 
tubes, but must not be taken as a basis for design throughout the whole range of 
dimensions employed in practice. 

The experiments of FAiUBAlRNjt on the other hand, were restricted to tubes of 
such relatively small length that he failed to realize the existence of a definite 
minimum below which the strength of a tube, however long, will not falL He also 
neglected the possibility of discontinuities in the curve of collapsing pressure at points 
where there is a change in the form of the distorted cross-section. In the light of 
these facts, figs. 3 and 4 help to explain his well-known formula, by which the 
collapsing pressure is given as inversely proportional to the length of the flue ; for a 
curve of hyperbolic form will represent as well as any other single curve the scattered 
points of fig. 4, and trial shows that the hyperbola 

is very closely an envelope of the discontinuous curve CBAE in fig. 3, in each case 
down to the point of least collapsing pressure. 

Validity of Investigation hy the Theory of Thin Shells, 

One important result of our investigation, which is apparently new, is shown by 
equation (83). It may be seen that collapse is practically dependent upon the pressure- 
difference alone, and that the absolute values of the pressures are immaterial. In 
view of this result, the objections raised by Basset against Bryan's treatment of the 
problem J require further consideration. 

These objections are : first, that the ordinary expressions for the stress- couples in a 
plate or shell, in terms of the curvature of its middle surface, are not valid when the 
surfaces are subject to pressure ; and secondly, that it is not legitimate to assume, as 
we must if sufficient equations are to be obtained, that the middle surface is 
unextended in a configuration of slight distortion. Hence the theory of thin shells 
is not applicable to this problem. 

The above difiiculties may be almost entirely overcome by a change in the method 
of investigation which is employed. It is customary to derive equations for the 
equilibrium of the distorted shell directly, and without reference to the position of 
equilibrium. Such procedure renders it necessary to make Bryan's assumptions, that 
the middle surface is unextended, and that the usual expressions for the stress- couples 

^ a/, footnote, p. 210. 
t Cf. footnote, p. 209, 
X See footnote, p. 210. 



ME. E. V. SOUTHWELL ON THE GENEEAL THEOEY OF ELASTIC STABILITY. 225 

are valid. But we may also proceed, as in the foregoing discussion, by first deter- 
mining the stress-system for the equilibrium position, and then deriving equations 
for an infinitesimal displacement. The stress-couples which appear in these equations 
will be due to the additional stresses introduced by the distortion, and since these, to 
a first approximation at least, vanish at the surfaces of the tube, they ivill be given 
with sufficient accuracy by the usual expressions. Moreover, when the distortion is 
two-dimensional (as in Basset's problem), the change in the '^ hoop " stress-resultant 
will be of an order which is negligible, so that the middle surface may be regarded as 
undergoing no extension relatively to the equilibrium position, even though its area 
may be sensibly changed in comparison with the unstrained configuration."^ 

The method of investigation just described, which follows the actual sequence of 
occurrences in the material, is suggested as in every way preferable to existing 
methods, for the investigation of any problem in elastic stability. For the present 
example, in particular, it leads to the same results as the more rigorous methods of 
this paper. 

Comparison with Existing Formidce, 

Previous discussion of the boiler-flue problem by analytical methods have, without 
exception, dealt with a tube subjected to pressure on one surface only, and almost all 
of them have been restricted to the case of an indefinitely long flue. Their results 
have, therefore, to be compared with our equation (80), when P2 is zero. It will be 
found that this equation agrees with the formula obtained by BryanI' and Basset :| 

Foppl's formula t omits the factor — , which measures the increased resistance to 

m^— 1 

flexure of a long tube as compared with a circular ring. 

The more general formula may be compared with that of Lorenz,J if ^2 be put 
equal to zero. It will be found that there is a serious want of agreement in regard 
to both terms in the expression (83). In support of the latter result, it may be 
urged that Lorenz' solution gives for the indefinitely long flue a result which does 
not agree with equation (80) (and, as we have just noticed, this is supported by 
previous investigations), and which vanishes, not when /c = 1, but when ^ = 0. Now 
the value 1, in the case of an infinitely long flue, corresponds to translation of the 
tube as a whole, without distortion, and the value to a change in the diameter 
of the tube, without any departure from circularity. It is clear that the applied 
pressures can have no tendency to maintain such a form of distortion, so that Lorenz' 
formula can hardly be correct. 

[■^ Added June 8. — The arguments of this section are more fully developed in a paper by the author 
"On the Collapse of Tubes by External Pressure," published in the 'Philosophical Magazine ' for May, 
1013 (pp. 687-698).] 

t Of, footnote, p. 209. 

X Cf, footnote, p. 210. 
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The ^^ Critical Length!' 

A, E. H. LovE'^' has investigated the rate at which the strengthening effect of 
circular ends falls off when the length of a boiler flue is increased. His result 
suggests that at ;a distance which is great compared with the quantity ^at the 
influence of the ends becomes negligible, and the flue collapses under sensibly the 
same pressure as a tube of infinite length ; hence, in order that '' collapse rings " may 
have any appreciable effect, their distance apart must not exceed some experimentally- 
determined multiple of this quantity. 

The greatest length of tube over which the ends exert any appreciable strengthening 
influence, or the least length for which collapse is possible under a pressure sensibly 
equal to the critical pressure, has been called by Prof. LovEf the '' critical length." 
It is a conception of great importance in experimental work; for, as we have seen, J 
tests on any length of tube in excess of this limit may be taken to give the strength 
of an infinite length of the same tube, and their results compared with the theoretical 
formula (80)§ : but as a basis for the spacing of '^collapse rings" it is superseded by 
the theory of this paper, which yields an expression for the greatest length of tube 
consistent with stabilit}^ when the thickness and diameter of the flue, and also the 
collapsing pressure, are given ; and Prof. Loye has suggested to the author that it 
would be better now to employ the term '^critical length" in this more general 
significance. As we have seen (p. 222), the length of the tube is some multiple of the 
quantity ajq, and we may therefore obtain from (83) the following formula : — • 

Critical length = ^ .- — — ^ . . (87) 



I 2E t ^w?-V >a} 

where M is a constant, depending upon the type of the collapse ring, and h has 
that integral value which gives the least value for the right-hand expression of 
equation (87). 

Before this subject is dismissed, it should be noticed that the theory of this paper 
does not support Prof. Love's estimate, mentioned above, of the rate of decay of end 
efiects. The term in equation (83) which depends upon the length of the tube may 
be regarded as negligible, compared with the constant term, when the ratio 



A 



2 /2 



* 'Proc. Lond. Math. Soc.,' XXIY. (1893), p. 208. 
t ' Theory of Elasticity ' (2nd edition), § 337 {h), 
X Page 224. 

§ In this sense the term " critical length" has also been employed by Carman, who began his research 
by investigating the strengthening effects of the end plugs with which he sealed his tubes for test. 
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has some sufficiently small value ; and -^ being inversely proportional to the length of 

the tube, we deduce for the '^ critical length," in the original sense of the term, an 
equation of the form 



l=/a/- 



3 
t 



where f is constant. Prof. Love, as has been said, has obtained an equation of the 

form 

L = /^ \/at, 

which is very different ; but he has informed the author that in the light of the above 
investigation (pp. 210-222) he does not regard his method as adequate/^ 

Solution for Tubular Strut: Speeial Case, 

We may obtain another simplification of the general determinant to ten rows by 
taking a zero value for h This corresponds to a type of distortion, possible in the 
case of a tubular strut, in which the axis remains straight and the cross-sections 
circular, the diameter varying in a sinusoidal manner. 

The ten-row determinant for this case is given on pp. 228 and 229 ; the factor —q^ 
has been cancelled from the sixth column, and terms in A have been omitted, so as to 
yield a result for tubes collapsed by end pressure alone. The expansion is only 
correct to terms of order r^, and for a first approximation we may also neglect the 
square and higher powers of B, which must be small in any case of practical 
importance. Investigating first the terms which are independent of r^ we obtain 

-i^ = ^^±lX (88) 

m q^ 

We may now employ the substitution 



m q' 



in the determinant, and expand it from the top row, neglecting terms of higher order 
than T^. 

A considerable amount of unnecessary labour may be avoided by a preliminary 
examination of the relative importance of the various terms involved. It will be 

P Added May J/.. — ^An argument in favour of the new formula may be drawn from physical considera- 
tions. The resistance offered by a tube to any given form of distortion is due partly to the extension 
and partly to the flexure which such distortion entails ; and it is clear that the relative importance of the 
extensional part increases as the thickness is reduced. Hence, other things being equal, the effects of the 
ends, which necessitate extension of the middle surface, are more important in a thin than in a thick 
tube; that is to say, they are sensible over a greater length.] 
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1 
found that B^ contains terms in (f and ™ , as well as terms independent of g. Thus 

the complete expression for B is of the form 

and it is clear that B has a minimum value when the axial wave-length has a finite 
value, given by 



2 4 r. m + 1 



yr'q^ =■ 2 — + ar . 

m 

This minimum value, which alone is of practical importance, is given, to a first 
approximation in terms of r, by the equation 



-^-'min. — ^'^ 'v/ — — ~°~~ y ^ V^/ 



so that the determination of a and ^ is not required. 
By expansion of the determinant we find 



r = 3 



2 "^ 



m— -1 ' 



and from (55) we deduce, for the minimum thrust required to produce collapse, 



<Smm. ~ SttEi^^A/ "3 . — (90) 

^ m — 1 



This expression is correct to terms in i\ 



Validity of Investigation by the Theory of Thin Shells. 

A complete investigation of the tubular strut problem must deal with lobed forms 
of deformation, since it is possible that one of these may require a smaller end~pressure 
for its maintenance than the circular form treated above. We have, therefore, to 
obtain a general expression for B (when A is zero) in terms both of k and q. 

The derivation of this expression, if we employ the rigorous methods of the present 
paper, will entail nothing less than the evaluation of the complete fifteen-row 
determinant ; for the existence of a ^' favourite type of distortion," of finite axial 
wave-length, which we have noticed in the particular case {k = O), is found by 
practical experiment to be equally a feature of the lobed forms of distortion, and 
shows that the terms in r^ are important. Now it will be shown that the value of 
Smin.) when k = 0, may be obtained, correctly to terms in t^, by the ordinary theory 
of thin shells ; and as there is no reason to believe that the latter theory will lead to 
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less accurate results when h has a finite value, it does not seem necessary to employ 
our more rigorous method, with the very laborious calculations which it entails. We 
shall therefore rely upon the approximate theory for the treatment of the tubular 
strut problem in its general form. Slight modifications in method will be introduced, 
as suggested above (pp. 224-225), and only the more important steps will be given 
here. 

Solution by the Theory of Thin Shells : General Case. 

We consider the stability of an element of the tube, originally bounded by the 

planes 

0, 6-f-§6, and z, z + Sz, 
as shown below — 

/ \ 

Perspective View ofEUmertt. 




Plan. View 
of ELemad.. 



Fig. 5. 

The other dimension of the element is the full thickness of the tube, denoted in 
this paper by 2t, The radius of the middle surface is a. 
The initial stress system is 

Pi = const. = — ■^~~ = [PJ (say). 

In the distorted position this system produces a radial force on the element, of 

amount 

1 



R 



[p J a se Sz, 



where E is the radius of curvature of a section of the distorted element by an axial 
plane (see fig. 5). 

It also produces a tangential force, in the direction of 6 increasing, of amount 



[PJ ai. se, 



where -i/^ (see fig. 5) = - — (l+e^^) Sz. 

a 00 
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The above system of distorting forces must be exactly balanced by the restoring 
system shown in the upper part of the figure. Hence we obtain the following 
equations of neutral stability : — 
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Now Pi and e^^ may be expressed in terms of the displacements of the middle 

surface, as follows : — 
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and the restoring system of stress resultants may also be expressed in terms of this 
system, as follows^' : — ■ 
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where D is the quantity 



m' 



m —1 



m^ . 



• • • 



V « 



(94) 



^ Cf, LovEj ' Mathematical Theory of Elasticity ' (second edition), Chap. XXIV. u', v\ w' have the 
same significance as in the earlier part of this paper, except that they now refer to the middle surface, and 
are functions only of 6^ and ^. 
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Eliminating Ti and Tg from equations (91), and substituting from (92) and (93), we 
have 
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where 
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Assuming a solution of the type (56), we find, as the criterion for neutral stability, 
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This equation, in its expanded form, is 
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Taking first the terms which are independent of t', and neglecting the square of ^ 
(which must be small), we find, as the first term in our solution, 
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When ^ = 0, this becomes 
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which shows that q must be great, if "^ has a value possible in practice. Similarly, 
when h>Q), we see that q must be small, and the approximate expression in this 
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We may now determine sufficiently approximate expressions for the terms in t^ja^, 
by treating q as great when ^ = 0, and as small when ^> 1. That is to say, we 
retain only the highest and the lowest powers of q in the two cases. "^ 
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When ^ > 1 , the important terms are 
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whence, to terms in t^/a" 
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with sufficient accuracy, when q is small. 
This leads to the result 
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(104) 



For practical purposes only the stationary values of S* are important. It is readily 
seen that the minimum value obtained from (l02) agrees with (90), and is therefore 



■^ In every case it is legitimate for practical purposes to neglect the term in ^^. 
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accurate as far as terms in f] we shall assume that (104) gives the same approxi- 
mation, which for practical purposes is quite sufficient. We then find, for values of h 
other than and 1, the expression 






(105) 



When ^ = 1, the axis does not remain straight after distortion of the tube has 
occurred. This is the type of distortion (sometimes called " primary flexure ") which 
was discussed by Euler, and it is easy to see that his result is identical with that of 
equation (104), which becomes in this case 



oV/g = 1 ' — Zir(Mfjli(^ , 



(106) 



The exact expression for the length of the tube, in terms of g, is not a matter of 
great importance in the present problem, because the wave-length corresponding to a 
minimum value of the collapsing pressure is in all cases small, and the strength of 
any strut of ordinary dimensions will therefore be given by equations (90) or (105), 
hito which the length does not enter. As in the case of the boiler-flue problem, we 
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Fig. 6. Strength of Tubular Struts. 

may illustrate the efiects of length upon the collapsing thrust by plotting the 
intensity of stress, or %/4:wat, against g"\ For this purpose we must take some 
definite value of the ratio t/a, and draw separate curves for difierent integral values 
of k The result is shown by fig. 6, in which the following values are assumed : — 
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From an inspection of these curves it is easily seen that as the axial wave-length 
increases the type of distortion which involves the least value for the collapsing 
thrust (and which the tube therefore tends naturally to assume) changes. For very 
short lengths we shall expect the circular type {k = O) ; then, as the length increasers, 
lobed forms of distortion, in which the value of k becomes less as the length increases. 
The limit is reached when k = 1 ; hence, the tendency of very long tubes is always to 
collapse in the manner discussed by Euler. 

It is also to be noticed that those parts of the different curves which lie to the 
right of their lowest points have no practical significance. The actual curve, which 
shows the effect of length uipon the value of the collapsing thrust, will approximate 
to the form shown in thick lines, since the wave-length (which varies as g~^) will 
naturally not increase beyond that value which involves the least collapsing thrust. 

Comparison with Existing For^mulw. 

The formulae of equations (90) and (104) may be compared with the results of other 
discussions of this problem. Equation (90) has been obtained by Lorenz,"^ and 

-— ^ is omitted, t 

The only existing solution for lobed forms of distortion is due to Lorenz,"^ and this 
is not in agreement with equation (104). In support of the latter formula it may be 
urged that Lorenz' formula does not agree with Euler's result when ^ = 1. 

It may also be remarked that the foregoing results for the tubular strut problem 
contradict Bryants theorem, that a closed shell cannot fail by instability, because 
distortion would involve extension of the middle surface ; for although the first 
terms in equations (l02) and (104) are due solely to extension of the middle 
surface, yet the compressive stress at collapse, as given by (90) or (105), may be 
insufficient to produce elastic breakdown in the position of equilibrium, if the ratio 

- has a sufficiently low value. 

a . . 

Stability of Tubes under Combined End and Surface Pressure, 

We shall not treat this case in any detail, but it requires notice in connection with 
the '^ localization of collapse '' which is observed in experiments conducted upon long 
tubes tested under hydrostatic pressure, the permanent distortion being generally 
confined to a portion only of the length of the tube. This result is not predicted by 
the theoretical formula (83), which suggests a steady fall in the value of the collapsing 
pressure as the wave-length increases ; and a partial explanation may possibly be 
found in the fact that the method of test has generally left a wholly or partially 

* Cf. footnote, p. 209. 

t For a similar omission in a solution of the boiler-flue problem cf. p. 225. 
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unbalanced end-thrust^ due to the water pressure acting upon the closed ends of the 
tube. 

It is clear that the expansion of the general fifteen-row determinant will give an 
equation of the form 

where a, ^, y... depend upon the dimensions of the tube and the type of the 

distortion. But in any practical case, as we have already observed, A and B must be 

very small quantities. It follows that an approximate solution may be obtained from 

the terms 

ot+^A + yB = (107) 

Let ^^ and <Sc be the values of the external pressure and of the end- thrust, each of 
which, acting alone, could produce collapse into the assumed type of distortion. Then 
equation (107) may clearly be written as follows : — 

fi=(l-f)ic, • . (108) 

where ^i and S are the values of the external pressure and end-thrust which can 
produce collapse when acting in conjunction. 

It may be seen from this equation that ^P^ can have a minimum value for some finite 
value of the axial wave-length when, and only when, & exists. If the end- thrust be 
entirely unbalanced, we have 

^= 7ra'(l-f-T)'fi, (109) 

and the collapsing pressure may, in this case, be determined from equation (108). 



Gbneeal Theory op Instability m Materials of Finite Strength. 

The Practical Value of a Theory of Instability. 

In the concluding section of this paper an attempt will be made to estimate the 
practical value of a theory of elastic instability ; to suggest ways in which we may 
hope to increase this value ; and to indicate the questions to which answers must be 
found in order that further advance may be possible. 

The first point which must be noticed is the non-realization in practice of our 
conception of a " critical loading,'^ owing to imperfections which always exist, and 
which violate our ideal assumptions. In any actual example the displacement of the 
system increases continuously with the load, and the system collapses at a smaller 
value of the load than our theory would dictate. It is necessary to inquire whether 
serious discrepancies are to be expected. 

In some mechanical problems the effects of imperfections may be calculated. We 
may take, as an example, the system illustrated in fig/1, and consider any one of the 
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many imperfections which occur in practice. For simplicity, let us assume that the 
sphere, bowl, and plunger are still smooth, rigid, and accurately formed, but that 
the line of thrust of the plunger is eccentric by an amount S, It is easy to see from 
fig. 7 that the displacement of the sphere from the line of thrust of the plunger, when 
the system is in equilibrium under a load P, is 



where 



d — r sin 6 — S + (R—7^) sin ^, 
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and these equations enable us to trace the steady increase of the sphere's displacement 
as the load on the plunger is increased from a zero value. 
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Thus in fig. 8 curves are drawn to connect P and d, for a value 3 of the ratio R/r, 
when the initial displacement ^ has the values 0, O'Olr, and O'lr respectively. At 
the points on these curves for which P has a maximum value, " collapse" will occur, 
since the equilibrium then becomes unstable. The locus of these points is shown in 
the figure by a broken line, and a dot-and-dash line shows the connection between ^ 
and the maximum value of P. From the latter curve it is evident that a small 
initial inaccuracy may cause a material reduction in the ^^ collapsing load"; never- 
theless the '' critical load" gives a limit which will be more and more nearly attained 
as our experimental accuracy is improved, and its investigation is by no means useless 
for practical purposes. 

When the problem is one of elastic stability, the discussion of imperfections by 
analytical methods will, in general, be beyond our power; but it is clear that similar 
remarks will apply. An '^ exchange of stabilities " at some " point of bifurcation "^ 
must be regarded as a purely ideal conception, and in practice there will always be a 
steady increase of distortion as the load is increased, owing principally to practical 
imperfections of form. A strut, for example, may be very accurately loaded, if 
suitable methods are employed, but its centre-line will never be quite straight ; the 
initial deflection which characterizes it may be regarded as composed of a series of 



■^ H. PoiNCARE, 'Acta Mathematica,' 7. (1885), p. 259. 
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harmonic terms, and when the load is applied one of these harmonics will be 
developed very much more than the others, just as one constituent harmonic may be 
developed by '' resonance '^ in an alternating current wave of irregular shape. In the 
ordinary strut problem this magnified harmonic is such that one-half wave occupies 
the length of the strut, but in other problems, such as that of the tubular strut, 
though there is always a '' favourite " or ^^ natural harmonic '' which is especially 
magnified, its relation to the dimensions may be more complicated."^ In any case the 
efiects of practical imperfections of form might be studied, if the analytical difficulties 
could be surmounted, by investigating the rate at which the amplitude of this 
^' natural harmonic '' increases with the load, when, its value in the initial configu- 
ration is given ; and the results of the investigation might be shown graphically by 
curves of distortion, similar in character to the curves of fig. 8, in which the abscissae 
represented the amplitude of the natural harmonic, and the ordinates represented the 
magnitude of the applied stress-system, or '' load." 

These '' curves of distortion" are of considerable utility for the study ot' problems 
in elastic stability, even though their true form can only be guessed. They help us, 
for example, to explain, and in some degree to remedy, the serious discrepancy 
existing between Eijler's theory and the results of experiments on short struts. 
The discrepancy has often been attributed to practical imperfections of form ; but it 
should hardly be necessary to point out that practical imperfections are likely to 
diminish rather than to increase in importance, as the dimensions of an elastic solid 
become more nearly comparable, so that they will never be more effective as causes of 
weakness than in struts of great length, which, as a matter of fact, give results in 
close agreement with Euler.'s formula. 

A more satisfactory explanation of this, and of similar discrepancies in other 
problems, may be found in the fact that the ordinary theory of elastic stability 
neglects the possibility of elastic break- down. If we attempt to draw '^ curves of 
distortion " for any single problem, we shall find that, apart from the other data of 
the problem, three possible cases exist, depending upon the elastic limit of the material 
under consideration : — 

(1) The material may be of infinite strength ; 

(2) Its elastic limit may be so high that the critical load, as determined by the 

theory of instability, is not sufficient to cause elastic break-down in the 
configuration of equilibrium ; 

(3) Elastic break-down may occur, even in the position of equilibrium, at a load less 

than the critical value. 

In the first case (which is, of course, purely ideal), the distortion due to loading 
will vanish when the loading is removed, and in this sense we may say that the 

* In the problem of the tubular strut, the " favourite harmonic " is, of course, defined by that value of q 
which corresponds to a minimum value of ^ in equations (102) or (104). 



240 MR. R. V. SOUTHWELL ON THE GENERAL THEORY OF ELASTIC STABILITY. 

material will never fail The '^ curves of distortion," if we could determine their true 
shape, would probably be approximately of the form shown in fig. 9. The theoretical 
methods of this paper enable us to fix the position of A, the '^ point of bifurcation," 
but give no information as to the form of AB, beyond the fact that it cuts OA at 
right angles.^ The other curves of the diagram will approach more and more closely 
the limiting form OAB as the initial value of the amplitude is decreased. 

In the second case, we have the additional complication of elastic break-down 
under finite stress, which reduces the resistance of the material and causes the new 
'^curves of distortion," shown by thick lines in fig. 10, to begin at certain points to 
fall away from the corresponding curves of fig. 9 (reproduced in fine lines for com- 
parison) ; these points will lie on some line such as CD, cutting OA at a point 
above A, and it is clear that to the right of CD the curves of distortion refer to 
displacements which do not wholly vanish when the load is removed. Total collapse 
of the system will obviously occur at the points of maximum load on the curves of 
distortion, and the locus of these points, which is shown on the diagram by the 
dot-and-dash line EF, may be termed the ''line of final collapse." 
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A knowledge of the true form of EF would enable us, when we are given the initial 
value of the amplitude, to predict the load at which the system will collapse ; and 
these quantities could be connected by another curve AG, which would show 
at once whether the resistance of the system to collapse is seriously reduced by 
practical inaccuracies of form. A complete theory of any problem in elastic stability 
must yield information on this very important point, as well as an expression for the 
" critical load " ; but in most cases more powerful methods would be needed for its 
derivation than are at present available. The investigation of the '' critical load" is 
therefore not without utility, for although never realized in practice, this forms a limit 
which should be fairly closely approached when considerable accuracy is possible. 

In our third case the " critical load," as deduced by theoretical methods, is more 
than sufficient to cause elastic break-down. We may proceed as before to draw 
hypothetical curves of distortion. The line G'D^ (fig. 11), which corresponds to the 

•^ It must not be assumed that AB is a horizontal straight line; in general, since the distorting 
effect of the applied stress-system, which varies as the deflection, increases less rapidly than the 
resistance, which varies as the curvature, AB will tend to rise from A. 
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line CD of fig. 10, will intersect OA at a point below A, and the other curves of 
distortion at correspondingly lower points. We have seen that the effect of local 
elastic breakdown upon fig. 10 was to deflect the curves of distortion from the forms 
which they would have assumed if the material had possessed indefinite strength ; 
and it is clear that this deflection will begin at lower values of the loads in the 
present case. We may therefore expect curves of the type shown in thick lines in 
fig. 11, where the curves already obtained are reproduced in fine lines for comparison. 
As before, we may draw a line A^F^ ''of final collapse" through the points of 
maximum load on the curves of distortion, and connect the collapsing load with the 
initial value of the amplitude by another curve A^G^ 

It is clear that the curves of distortion must tend to a limit which is no longer 
OAH, but some other curve OA^H^, where OA^, the critical load under the new 
conditions, is more than sufficient to produce elastic breakdown, but less than OA. 
We can see further that the curve A^G^ is not likely to fall away from A^ much more 
steeply than AG from A in fig. 10. The great weakness of short struts in practice, 
compared with Euler's theoretical estimate, is now explained. Whereas long struts 
come within the conditions of fig. 10, the failure of short struts will be repre- 
sented by fig. 11, and occurs at comparatively low stresses, not because practical 
imperfections have a greater effect upon the strength, but because OA^, the true 
value of the critical load, is less than OA, the value which Euler's theory would 
dictate.^ 

It is the rule, rather than the exception, that the critical load, as found by the 
ordinary theory of elastic stability, is more than sufficient in practice to produce 
elastic break-down. This may be readily seen in reference to any particular example. 
In the case of the tubular strut, fig. 10 is only applicable when the ratio of diameter to 
thickness is greater than 560 (for an average quality of mild steel), and for thicker 
tubes the critical load falls, apparently by a very considerable amount, t below the 
theoretical estimate. The determination of the critical load, in cases where this is 
more than sufficient to produce elastic break-down, is thus a problem of great 
importance, since it forms a limit which can never, under any circumstances, be 
exceeded. In the ordinary strut problem the determination can be effected without 
difficulty, and an apparently new field is thus indicated for research. The distin- 
guishing feature of its problems is the dependence of the stress-strain relations upon 
the past history of the material, rendering absolutely necessary a method which 
follows the actual cycle of events up to the occurrence of collapse. 

[* Added May 11. — Since this paper was written, the author's attention has been drawn to a 
dissertation by T. voN Karman ('Untersuchungen iiber Knic'kfestigkeit/ Berlin, 1909), in which the 
forms of these " curves of distortion," for solid struts of practical dimensions, are deduced both from theory 
and from experiments. Karman also gives a relation equivalent to that of equation (112).] 

t Experiments conducted by the author upon seamless steel tubes showed failure under loads which 
were in every case little more than sufficient to produce " permanent set." 

VOL. COXIII. — A. 2 I 
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Stability of Short Struts. 

This problem has been discussed elsewhere by the author,^ and it will be noticed 
here only at sufficient length to indicate the directions in which further research is 
needed. We have to derive an expression for the collapsing load of a straight struts 
when this is more than sufHcient to cause elastic break-down of the material ; and we 
proceed as before by considering three configurations of the strut : (l) before strain ; 
(2) in a position of neutral equilibrium under uniform end- thrust ; and (3) in a 
position of infinitesimal distortion from the second configuration. 

For a first approximation we may say that cross-sections remain plane in the third 
configuration, so that the diagram of longitudinal compressive strain for any cross- 
section is as shown in the upper part of fig. 12 ; the horizontal line fg shows the 
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uniform strain of the second confi.guration. Then, if fig. 13 be the stress-strain 
diagram for a compression test of our material, and this uniform strain corresponds to 
a stress p which is represented by the point B, we see that to the right of the point 
F in fig. 12 the longitudinal compressive stress in the third configuration must be 
greater, and to the left less than p. 

Now it is a well-known property of metals that if at any point B on the stress- 
strain diagram, beyond the elastic limit, we begin to decrease the load, the diagram 
is not retraced, but that we obtain a line BC which is parallel to OA.t It follows 
that the ratio decrease of stress 

decrease of strain 

is still given by E, Young's Modulus for the material. On the other hand, the 
diagram shows that if we increase the load beyond B by an infinitesimal amount, the 

"^ • increase of stress 

increase of strain 

is a smaller quantity E^, which may be found from the slope of the diagram at B. 

* ' Engineering,' August 23, 1912. 

t A. MoRLEY, * Strength of Materials,' § 42. 
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We are considering an infinitesimal distortion in the third position, so that if we 
represent the increase of strain in fig, 12 by \z^ the increase of longitudinal com- 
pressive stress to the right of F may be taken as WXz, and the decrease of this stress 
to the left of F as EX^;. Hence we obtain, for the section under consideration, the 
diagramx of longitudinal stress which is shown in the lower part of fig. 12. The 
uniform stress of the second, configuration is shown by the horizontal line In^ and it is 
a condition for neutral stability in the second configuration that no increase of thrust 
shall be required to maintain the distortion. If the cross-section of the strut is 
rectangular, of dimensions a x 2^, it follows that the triangles Imh and qmn must be 
equal in area, or 

This relation fixes the position of F on the cross-section, and in terms of the 
dimensions shown in fig. 12 we may write for the moment of resistance about G, 

'¥A^a\z{z—h)dz-\- 'EiaXz{z—h)dz 

= faXE^^+iaX (E--EO {¥-^^-21% 

But if y is the deflection of the strut at the point G, in the infinitesimal 
distortion, we have, as in the ordinary theory of bending, 

where x denotes the distance of the section from one end ; and for equilibrium in the 
third configuration M must be equal to the bending moment due to thrust, or 2atpy : 
hence, 

i>t/+gxiE^^|l+i(l-|)(|-.3^-2)| = (112) 

Equation (112) shows the modification which must be introduced, to take account 
of elastic break-down, into Eulek's equation 

py+~^x^'Eif =^ 0, . (113) 

and it is easy to see that if Z is the length, calculated from (113), of a strut which 
can just support the stress p, and V the length as calculated from (112), then 



But, from (111)5 
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so that, finally, 

— - ( 1 1 a\ 

— *— . "v~^ • * • • * • • * * • \X.i jt / 

1+ /y ^ 



This result leads to a simple method by which the collapsing loads of short struts 
may be obtained graphically from the compressive stress-strain diagram. A full 
explanation is given in the paper to which reference has already been made, and a 
comparison, showing satisfactory agreement, is made with the results of experiments. 
One conclusion of some practical importance may be noticed : the curves of collapsing 
stress show that great ductility of material is by no means desirable in struts, 
the primary requisite being a high elastic limit. 

Need for Further Research, Conclusion, 

With slight modification, the theory just given for short struts might be applied to 
the problem of circular rings under radial pressure ; but these appear to be the only 
cases in which we can at present discuss the stability of overstrained material. In 
any problem dealing with plates or shells distortion from the equilibrium position 
must introduce new stresses, in directions perpendicular to that of the stress which 
has caused elastic failure. The circular type of distortion in a tubular strut, for 
example, will introduce " hoop " stresses, and at present we have no knowledge of the 
corresponding stress-strain relations when "set" has occurred. 

This and many other stability problems may be regarded as special cases of a 
general problem, viz., the determination of the changes of strain which occur when an 
infinitesimal stress-system, defined by principal stresses g, r, 5, is impressed upon a 
material already overstrained by a simple stress p. The problem is not simple, and 
its solution would probably entail much theoretical and experimental work ; but this 
would be justified by the importance, both for theory and practice, of its applications. 

In conclusion, the author desires to express his indebtedness to Profs. Love and 
HoPKiNSON, for valuable criticism and advice ; to Mr. L, S. Palmee, for the photo- 
graphs reproduced in fig. 2 ; and to Messrs. H. J. Howabd and D. P. Scott, for 
assistance in the prosecution of the experiments described on p. 223. He also takes 
this opportunity of thanking Messrs. Stewarts and Lloyds, Ltd., of Glasgow, for 
gifts of very accurate steel tube for experimental purposes. 
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